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Abstract 

The relation between the correlation energy and the entanglement 
is analytically constructed for the Moshinsky's model of two coupled 
harmonic oscillators. It turns out that the two quantities are far to he 
proportional, even at very small couplings. A comparison is made also 
with the 2-point Ising model. 



1 Introduction 

The concept of entanglement has been recently considered by many authors 
[1] in connection with several properties of the quantum systems and as 
a potential resource in quantum computation and information processing. 
Moreover, entanglement has also been recognized to play an important role 
in the study of many particle quantum systems [2] and experimental mea- 
surement have demonstrated that it can affect macroscopic properties of the 
condensed matter [3]. However, a new area of research has been opened by 
[4], where it was argued that the entanglement (a quantum observable) can 
be used in evaluating the so-called correlation energy: that is the difference 
between the true eigenvalue energy of a given molecular system, composed 
by identical components, with respect to that one prescribed by the Hartree 
- Fock (HF) approximation method. In [5] the case of the formation of 
the Hydrogen molecule was discussed, and a qualitative agreement between 
the von Neumann entropy of either atom (measuring the entanglement of 
formation of the whole system) and the correlation energy as functions of 
the inter-nuclear distance was shown. However, the extension of this idea 
to multi-atomic molecules and its effectiveness remains still unclear [6,7]. 
The correlation energy in this context is an artifact of the approximation 
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procedure, then it is not a physical observable and, by second, it can be 
modified by the adopted method of calculations. In this sense the conjec- 
ture of [5] is not a well posed question and the answer to it must be only 
qualitative. Nevertheless, folk says that still entanglement has to play some 
role in taking into account the difference between the approximated com- 
putation of the (fundamental) energy level for a typical factorized wave 
function, as provided by the Hartree-Fock approximation, and the "true" 
energy eigenvalue. 

This paper is devoted to a deeper analysis of it, trying to find an oper- 
ative and quantitative estimation of the entanglement effects on the energy 
level evaluation by approximated methods. Better, one would like to clarify 
if the correlation energy is in a one-to-one correspondence with the entangle- 
ment of the fundamental state of bipartite systems. To this aim we analyze a 
very simple model of coupled harmonic oscillators introduced by M. Moshin- 
sky about 40 years ago in the paper [8], and try to arrive to a quantitative 
measurement of the entanglement effects on the correlation energy. 

In Sec. 2 we briefly review the Moshinsky's model, pointing out the 
goodness of the HF approximation for the fundamental state, both con- 
sidering the correlation energy and the fidelity with respect to the exact 
solution and for weak couplings of the two oscillators. In Sec. 3 we evaluate 
the von Neumann entropy, which is the unique measure of entanglement for 
pure states, for the fundamental state of the Moshinsky's model, by trac- 
ing out the density matrix operator w.r.t. one of the component oscillator. 
Eigenvalues of the reduced density matrix are computed by discussing cer- 
tain integral equations. In order to compare the von Neumann entropy and 
the correlation energy, in Sec. 4 we consider the minimized mean squared 
deviations of the two quantities. A discussion of the relative deviations is 
also performed. Furthermore an analytic expression of the entropy in terms 
of the correlation energy is provided. A comparison with the analogous ex- 
pressions for the two point Ising model is discussed. Asymptotic formulae 
for very small couplings ( or correlation energy) are considered at the end. 
Some final remarks are addressed in the Conclusions. 

2 The Moshinsky's Model 

In order to evaluate the efficiency and how good the HF mean field method 
is in computing quantum states, Moshinsky in [8] proposed a very simple, 
but non trivial, model of two coupled spin-^ harmonic oscillators in 3 di- 
mensions, in which all computations can be analytically performed up to 
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the end. 

In adimensional unities, the Hamiltonian of the model reads 

1 (^2 ^2 -^2 ^2\ 1 ^ \ 2 

H = ^[pi + P2+ri+r2)+-Kl^ri-r2) , (1) 

where and denote the position and the momentum operators of the 
i-th particle, respectively. The constant K parametrizes the interaction 
strength of a supplementary harmonic coupling between the two oscillators. 
One can interprete such as a model for two identical interacting atoms (or 
molecules). Thus, the coupling woTild represent the interatomic interac- 
tion, which one can conjecture weaker than the electron-nucleus interaction. 
Thus, we will dwell upon the case < K < 1 and the coupling term in 
(1) will be considered as a perturbative term. So changing the value of K 
may be interpreted as a change into the interatomic interaction. Finally, 
let us observe that since we are interested into the fundamental state, the 
spinorial aspect of the problem is not relevant. In fact, the total spin state 
must be singlet. For self-consistency of the present paper, here review the 
Moshinsky's results. 

The energy spectrum of the system can be easily computed in the form 



En,m = ^ (l + (1 + 2X)V2) +n + m{l + 2X)V2, (2) 

where m and n are positive integers or 0. For any K the spectrum shows 
degeneracy, partially broken for K > 0. However, the lowest level is 
always simple. Moreover, eigenvalue crossings occur for higher eigenvalues 

at isolated points of K . 

The normalized wave function of the fundamental level is given by 
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where the mean and relative position have been defined, R = ^ ^ and r 

= ^ ^ respectively. Looking at the spin state of the system one can see 
that it is always entangled. So the question is if also the positional factor 
has the properties of an entangled state. 

Applying the standard HF mean field approximation to the ground state, 
one is led to the wave function 
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Figure 1: Exact and HF approximated fundamental level (on the left). Cor- 
relation Energy (on the right) 



^HF r) ^ = ,-3/2(1 +^)3/4^-(l+i.)V.(«.+..)/2 I 1, l,0,o) , (4) 

corresponding to the approximated eigenvalue 

EnF = ^{l + K)^'^. (5) 
Defining the correlation energy as 
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^corr = ^HF " ^0,0 = + K + Vl + 2Kj , (6) 

one obtains a first indication of how good the HF approximation is. For 
simplicity this is reported in Fig. [TJ Of course, by Ritz's theorem one has 
-E'hf > -£^0,0 5 but the approximation looks very good for small K, indeed 
they differ at the second order near iT = as 

Ecorr=^-^-^-^+0{K'), (7) 

SO that the relative error is less than 3.5% in the range of < -fC < 1. 
Moreover, the explicit expression of the overlap ( or the squared fidelity) of 
the exact and the HF wave function |(^'hf | ^o)P as a function of K is 

{K + {i + vttk) (1 + VI + 2K)y 
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the graph of which is shown in Fig. [2j It is immediate to see that the 
distance between the two wave-functions is less than 95% for < K < 1 . 
Thus, one can figure out that adding to the HF state further corrections 
surely improve the estimations of the energy eigenvalue and increase the 
fidelity, but the simplest factorized expression in (4) will be lost. Indeed, 

this is the case if one considers the exact expression of '^oiR,^) in (3)- 
Qualitatively this means that the two oscillators in the correct fundamental 
state are entangled, contrary to what happens in the approximated state. 
Thus, one would like to attribute the mismatch between the energy of the 
two states (i.e. i^corr) to the entanglement properties. We will analyze this 
idea in the next Section. 



3 Entanglement Estimation 

As it is well known, the main estimator of the entanglement for pure states 
is the von Neumann entropy of the reduced density matrix pi of one of the 
component subsystems, i.e. 

5'vN [pi] = -Tr [pi log2 pi] , with pi = Tr2 [p] (9) 

where in the second relation one has traced out only with respect to the 
second subsystem, when a proper basis of states has been chosen in the 
Hilbert space of the whole system. In the position - spin representation the 
complete system density matrix p for the fundamental state (3) is given by 
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p{ ri, r2,r i,r 2 , (10) 



P ( r 1, r 2, Ay 2) = TT-^l + 2K)3/4e-(«^+«'^)/2e-(i+2i^)^/^ (.^+r.-)/2^ 

(11) 

where one has introduced the supplementary variables R' = ^ ^ and r' 



An analogous expression holds for the |*hf) state (4). The main dif- 
ference between the two functions is that in the latter case the density 
matrix reduces in just the product of gaussian normal distributions in the 

real variables ( r i,r 2,r' i,r'2), with the same variance. Instead, in the ease 
of density matrix (10) the different coefficients in front of {E^+R'"^) and 
(r^+r'^) make the presence of mixed (entangling) terms of the form ri-r2 

and r'i-r'2- This has consequences also at level of one particle space 
distribution density, which is given by 



Poi (r ] = p(r,r2,r,r2]dr2 = v^k+i+i ^ 

^ ^ J ^ ) 7r3/2(V2^TT+l) ^ 

(12) 

which differs significantly from the one particle density obtained from the 
HF state, i.e. 



_ (K + 1)3/4 
PHF7 = e ^ ^ , {V6) 

as, for instance, it can be graphically seen in Fig. 3 for a particular value of 
K. Thus, in the exact state the mean squared deviation in the position mea- 
surements of one particle, regardless the other one, results greater than in 

1 /2 

the approximated case by a factor ^ ^^"'^^^1+2^^^^'^ "^ ' Correspondingly, 
the mean squared deviations for the modulus of the linear momentum are 
smaller, by the inverse of the above factor, in the exact case with respect to 
the HF one. 
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Figure 3: Exact and approximated one particle probability densities, for 
K = 1. 



The analytic 1-particle spacial density matrix operator can be given in 
position representation pi =/ pi , r'^ • dr' hy the kernel 

^\ _ (23/^(2i^+l)3/4) 



exp 



(a/ZKTT-I) ^•r'-(j!:+3V2K+l+l) ( + r' 
4(^2^+1+1) 



(14) 



which we have to use into expression (9) in order to compute the corre- 
sponding von Neumann Entropy. Actually, one has to take account also of 
the spin state. But, a well known property of the von Neumann Entropy 
says that 

5'vN [/5 (t] = 5vN [p] + 5'vN [^] , (15) 

for any factorized density operator p a . Since this is the case for the 
system we are considering, where the reduced matrix in the spin space state 
is ^ 1 proportional to the identity matrix, the von Neumann Entropy for 
the one-particle subsystem can be defined modulo a constant term equal to 
1. 

Furthermore, the expression (9) becomes S'vNiPi] = —J2iPi^'^92 (Pi) 
terms of the eigenvalues /Xj of the operator pi, which is acting on the infi- 
nite dimensional space of squared sommable functions. So, the eigenvalue 
problem for pi can be written in the form 
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J Pi ^ r , r'^ u (^^''^ = /X u ^ r ^ . (16) 

As one can see from Eq. (14), the kernel of this integral equation is 
symmetric. Moreover, it is of Hilbert-Schmidt type, since the coefficient of 

r' is negative. So the spectrum is real and discrete. Finally, in cartesian 
coordinates it can be decomposed in the product of three one-dimensional 
kernels with the same gaussian structure 



/4 



(v^27?+T+l)''"7ri 



1/2 



exp 



4(^2^+1+1) 



(17) 



Accordingly, the eigenf unctions of Eq. (16) can be factorized in the product 
of three functions, each of them depending only on one real variable. Thus 
we put for the eigenfunctions and the eigenvalues 

Ul,m,n (^r^ = w{x)w„,{y)Wn{z), /i/,m,„ = t'/t'mJ^n, (18) 

SO that the problem is reduced to solve the integral spectral problem in one 
dimension 



P(X,x')wi{x')dx' = vi wi{x). (19) 
Recalling the standard formulae for gaussian integrals [9] 



/oo 
-oo 



xV-'--'^dx' = -l-y|9-(ge1^), (20) 



it is natural to look for solutions of Eq. (19) of the form 

wi{x)=Pi{x)e-^'^, (21) 

where P/(x) is a polynomial of degree / in the variable x and 5 is a real 
positive coefficient. 

Substitution of the expression (21) in Eq. (19) leads to a pure polynomial 
expression iff one sets 

6=^{l + 2Ky/^. (22) 
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Figure 4: The eigenvalues fiQ and /xi of the reduced density operator pi. 



Then, Eq. (19) reduces to (^Z+iJ-dimensional hnear eigenvalue problem 
for the coefficients of x in the polynomial Pi{x). However, the particular 
form of the coefficients in the r.h.s. of formula (20) makes this problem 
triangular, so one can easily provide the non degenerate spectrum 



where 

C 



yi = Cc\ (23) 
2x/2(l + 2Kfl^ 



l + K + 2{l + 2K) 1/4 + 3V1 + 2K + 2(1 + 2i^)3/4 

l + K -y/l + 2K 

^ ~ l + K + 2{l + 2i^) 1/4 + 3V1 + 2K + 2(1 + 2K)3/4 ' 

Of course, these eigenvalues are positive and their sum is equal to 1, be- 
cause related to a matrix density operator. On the other hand, accordingly 
to (18) the eigenvalues of the one-particle density matrix (14) are given by 

= C^c '+"^+" = C^c^ e N U {0}, (24) 
thus the eigenvalue /x^ is degenerate of order 

deg[/Xfc] = y + y + 1. (25) 

For comparison in Fig. 4 we report the graphs of a-^d /ii. Hence, if 
we are allowed to interprete the eigenvalues of the density matrix operator 
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Figure 5: The traee of pf, as a function of the couphng constat K. The 
state of one the Moshinsky's oscillators becomes more mixed for increasing 
K. 



pi as the probabilities to find the one particle subsystem in the states of a 
iT-parametrized family of harmonic oscillators, for the fundamental one it 
is very close to 1 and slowly decreasing in K. At K 1, the probability to 
find it in the first excited state is about 1.76% and exponentially smaller for 
the highest excited states. The lack of coherence can be estimated also by 
computing the Trace [pf]. Indeed this quantity is 1 only if it corresponds to 
a pure state, but in the present case one has 

, _ 8(1 + 2Kf/^ 

{l + VTT2Kf ^''^ 

which is a decreasing function in K, as can be seen from Fig. 5. 

Now, we are in the condition to compute the von Neumann entropy, by 
rewriting formula (9) as 

'S'vN [pi] = - J2T=o deg [Pk] Pk log2 Pk = 

(27) 

{2 (1 + x') In [1 + x] - 2x In [4x] - (-1 + xf In [-1 + x'] } X"', 



3 

top] 
where 



X = (1 + 2K)4. (28) 
A graphic of -SvNipi] is given in Fig. 6 for the range of < K < 1. 
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Although in the present case the entropy is not upper bounded, since 
we arc deaUng with a system with an infinite-dimensional Hilbert space, 
the expression (27) is very to the analogous one for the 2-point Ising model 
[5,10], for which the 1-particle von Neumann entropy and the correlation 
energy are 



rt2— Ising 



[P 



{ (-2 + V4 + A^) Zn 



1 



+ f 2 + \/4 + A^) In 



(29) 



£^corr = V4 + A2 - 2 , 



2 + ViW 



(30) 



where A is the normalized ferromagnetic coupling. 

For large K the entropy SyN [pi] increases logarithmically according to 
the expansion at infinity 



•SvN [Pl 



?,ln[K] 
4/n[2] 



21 3 
T ^ ln{2) 



+ 



1 



On the other hand, the behavior of SyN [pi] near K = can be described 
by its series expansion 



•SvN [Pl 



3 3 

4 ^ 16/n(2) 



3^3 



+ 



177 



189 



512 /n(2) 



3ii:2 (16 - 32K + 59^2) 
128Zn(2) 



Zn(i^) + O [K^] , (31) 
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but this approximation becomes inaccurate very rapidly. Prom the above 

expression one sees the asymptotic behavior of the entropy close to is 
controlled by a logarithmic term, differently from the correlation energy 
(6), which has a pure power expansion. Then, we cannot expect a great 
similarity between the two functions, also at very small values of K. 



4 Comparing Entropy and Correlation Energy 

The comparison of entropy and correlation energy requires at least a common 
scale of unities, the former being a pure number. So, the simplest proposal 
is to study a family of deviation functions of the form 

A{a,K) = aS^n[pi]-Ecorr, (32) 

where a is a parameter to be properly chosen. The choice of a can be done 
in several ways but, trying to minimize the values of A (a) in a range of 
K, a criterion could be given by the minimization of the squared deviations 
functional 



I[a] = [ A{a,K)^dK. 
Jo 



(33) 



Quite arbitrarily, the domain of integration is determined by our interest to 
small perturbations of the 2-independent harmonic oscillators system. This 
is equivalent to set 

«min = .T" ^ 0.318949, (34) 

Jo -^^vN [Pij dK 

where a numerical evaluations is needed because of the intricate form of the 
integrand functions. The graphic of the function A(ainin,K) is given in Fig. 
7 . An estimation of the relative deviations of the correlation energy with 
respect to the entropy is given in Fig. 7. In both cases we have considered 
the relative deviation at the minimum of the squared deviation functional 
(31). As one can see, in both cases the relative error is quite big, particularly 
close to 0, where the behavior of these functions are largely controlled by 
logarithmic or cuspidal singularities. In particular, while in the expansion 
of A(Q;niiinK)/ SyTSf [pi] a finite cuspidal singularity appears at 0, i.e. 

|A(a.m,K)| _o.o35439+o(i^), A f \M-r..r., K)\\ -8.20851 x 10- 



SyN [Pi] ■ ' ' dK \ Sy^ [pi] J 
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A(ami„,K) 




Figure 7: A plot of the minimized deviations between entropy and correla- 
tion energy. 



|A(a„in,K) |/SvN [^oi] 




Figure 8: The relative minimized deviations between entropy and correlation 
energy. 
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while the relative deviation with respect to the correlation energy shows a 
true logarithmic divergence as 



|A (a^in, ^ 1 + 114945 ^ ^g-ie ^ ^^^^^ ^35) 



corr 



making very questionable the meaning of these quantity, taking into account 
also the tiny absolute values of the expansion coefficients. 

On the other hand, correlation energy and entropy look almost pro- 
portional in the restricted range ~ 0.6< K < 0.8. In this region of the 
K parameter it is confirmed in explicit way how the correlation energy is 
a direct effect of the entanglement. In general such a relation is less obvi- 
ous. In fact, first let us observe that both functions S'vN [pi] and i^corr are 
both monotonically increasing in K. Therefore, since .Ecorr is a quite simple 
algebraic function on the coupling constant, it can be easily inverted and 
replaced into .Svn [pi] , obtaining again a one-to-one correspondence, say 
SvN (-E'corr), the cxpression of which is a little involved, actually 

-SvN (^^corr) = ^^=]^ {4V3(r + 3) In ( V6f+9 + 3) (36) 

-2V3ln{6) (r + 3) - 2 (V3(r + 3) - 3^2 r + 3) /n(r) 
-3V2r + 3 Zn(2r + 3) - 3V2r + 3 /n(12)} , 

where we have introduced the auxiliary variable r = Ecorr+V^^/ Ecorr {^corr + 3). 
A plot of this function is given in Fig. 9. The analogous expression for the 
2-point Ising model results directly by substituting Ecorr given in Eq. (30) 
into the expression (29) for the entanglement. A plot of this function is as 
is also sketched in Fig. 9. 

In both cases the entanglement is an increasing function of the correlation 
energy. In particular, one can search explicit but approximated expressions 
of the entanglement for small values of Ecorr, corresponding to small cou- 
plings. Indeed, including logarithmic corrections at the lowest order near 
-Ecorr ~ 0, one obtains 



5 (J, {l+ln[6]-ln[E,,,r]) ^ 

'-'vN (.-C/corr j ~ 27n[2] -^corrj 

$2-Ising . p X {l + 2ln[2]-ln [.Ecorr]) ^ 

{^coii) ~ 4,ln[2] -f^corr, [OO) 

for the Moshinsky's oscillators and the 2-point Ising model, respectively. As 
one can see, the two asymptotic formulae above are very similar, but not 
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Figure 9: The entanglement as a function of the correlation energy for 
the Moshinsky's model (continuous curve) and for the 2-point Ising model 
(dashed curve). 



sufficient to argue that there exists a general argument allowing to compute 
the coefficients appearing in the above developments. 

For large perturbations, or equivalently for -Ecorr oo, one has the 
asymptotic expansions 

S,. (E,„,) ^ ^ + -"-■W - '"'"'> . . O (-^) (3.) 

ln{4) 2ln(2) \Ecorr J 

for the Moshinsky's model, while for the finite levels 2-point Ising model the 
entanglement will approach the limiting value 1 as 



-'corr ' 

Concerning the singular behavior of the entanglement as a function of 
the correlation energy it does not seem related to the specific way of its 
estimation. In fact, in alternative to the von Neumann entropy one can 
use the concurrence [10], which uniquely defines the entropy at least for the 
Ising model (see Eq. [29|) . In this case the concurrence takes the form 

^2-/sm<; ^ ^ (41) 



Furthermore, expressed in terms of the correlation energy in (j30p . it becomes 



^2-Ising _ \/-5'corr (-E'corr + 4) ^^^^ 
-^corr ~t~ 2 
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As one can see, the first derivative of this expression contains an algebraic 
singularity at instead of a logarithmic one, as it shown by Eq. (|42p . and 
it has a monotonic algebraic increasing to the limiting value 1. 

Inspired by [10], one can define the analogous of the concurrence for 
the Moshinsky's system (at least limiting ourselves to the set of considered 
ground states) as 

C=1-Tt[pII (43) 

which takes values in the range [0, 1[ . It is invariant under local unitary 
transformations on the separate oscillators (reduced to changes of phases) 
and provides a unique mapping for 5vN [pi], by using Eq. (28) into (26) 
and, then, replacing into (27). 

As a function of the correlation energy, the above defined quantity takes 
the form 

/ \3/2 

3\/3 I lEcQ-pf + lyl^ E(^orr {,Ecorr + 3) + 3 J 

C{E,orr) = 1 ) ^3^, (44) 

Ecorr ~\~ \/2a^/ Ecorr {-^corr ~l~ 3) -|- 3 



which is regular in the origin, but it is not its second derivative. Again a 
singularity is signaling the a greater rate of increase of the entanglement 
with respect to the correlation energy, for small values of the latter. 

One advantage of the expression ([441) is that it can be easily inverted, 
providing a quite simple expression of the correlation energy of the Moshin- 
sky's model in terms of the entanglement, which was one of the original 
motivations of the present work. Specifically, one has the correlation energy 
as a function of the concurrence (41) 



-i^corr — (^i_Q^5/6 



2C - ^/r^ 2^(1 - C)4/3 - 3(1 - C)2/3 + 2 - 3 -2 



+ (1 - C)i/6 ^(1 _ C)4/3 _ 3(1 _ C)2/3 + 2-1 . (45) 



In particular, the power expansion of the above expression near concur- 
rence is given by 

Ecorr - C + C'/' + ^ + (46) 

This expression gives a direct relation between the correlation energy and 
the entanglement via the quantity defined in (j43p . much simpler than by 
using the von Neumann entropy. 
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Figure 10: The overlap for HF w.r.t. the exact ground state of the Moshin- 
sky's model entanglement as a function of the 



Finally, the fidelity of the fundamental state of the Moshinsky's model 
with the corresponding HF state, or equivalently the overlap (8), can ex- 
pressed as a function of the entanglement. In some sense, we are comparing 
two different way to measure the "distance" between the two curves of states, 
even if neither quantities actually have the properties of a distance. How- 
ever, also in this case a monotonic function can be obtained for any pair of 
states corresponding to the same coupling constant K, or correlation energy 
Ecorr- But, in general, an explicit expression of such a function is unknown 
for generic states of the Moshinsky's model. A plot of the function is shown 
in Fig. [TOl Also for this particular situation, even if regular, the behavior 
at the origin possesses a singularity in its higher derivatives. 



5 Conclusions 

In the present paper, from point of view of the entanglement theory of quan- 
tum systems, we have described the relations occurring between two special 
curves in the space states of a family of bipartite system (the Moshinsky's 
model of harmonic oscillators), continuously parametrized by a coupling 
constant. To varying of that, one curve contains the exact ground states, 
the other the states provided by the HF method approximated solutions. Of 
course, for vanishing coupling the two curves flow out from the same state, 
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but we would characterize how they separate in terms of entanglement. The 
peculiarity of the latter curve is to lie always in the manifold of entan- 
glement, while the other one suffers a monotonic increasing entanglement. 
This could provide a measure of the "distance" between the two curves and, 
maybe, of the geometry of the space state around them. On the other hand, 
a similar, but not necessarily equivalent description is given in terms of 
the correlation energy, which in principle is defined only for pairs of corre- 
sponding states (at same coupling K) in the two curves. We have proved 
that entanglement and correlation energy are one-to-one along these curves. 
However some peculiarities arise. First, they are far to be proportional and 
only for certain intervals of the coupling constant their squared deviations 
can be considered small to some percentage. Second, at origin they have 
a quite different rate of increasing. We show that this phenomenon occurs 
not only in terms of the von Neumann entropy, but also by introducing an 
adaptation of the concurrence. From the view point of physico-chemical 
calculations the above observations say that an artifact of the calculation 
methods, the correlation energy, can be interpreted and calculated in terms 
of the entanglement of the wave function of the exact solution. Moreover, 
we have quite simple algebraic relations in terms of the concurrence (see Eq. 
()46p ). However, at the moment we have not a general method to compute 
directly the coefficients of such type of expansions. These could be very 
useful in order to have an alternative a priori estimation of the error made 
in computing the correct expectation values of the energy. 

We have shown also that the overlap of the exact wave functions with the 
corresponding HF approximations are in a one-to-one correspondence with 
the correlation energy and the entanglement. In particular a monotonic 
decreasing of the overlap occurs as a function of the entanglement, with 
a non analytic behavior at the origin. Differently as above, such a type of 
relation may be useful in the estimation of the entanglement, which could be 
quite complicated to compute for bipartite systems with high inner degrees 
of freedom, as the dimers of complex molecules (see [6-7] for instance). 

Finally, we would stress that if a general principle of nature would state 
that it is impossible to create (or increase) entanglement between remote 
quantum systems by local operations [11] (in analogy with the second law 
of the thermodynamics), expressions like (32) may suggest proper modifica- 
tions in presence of globally controlled operations, analogous to the isother- 
mal transformations. In this respect, the true meaning of the minimizing 
parameter a is still obscure and could be a further direction of investigation. 
It could be just a case but, as noticed in [7], the properties of the function 
A {aram,K) could be used to predict special configurations of the considered 
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physical systems, in a sort of balance between entanglement and energy. 
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